DESY 99-166 
November 1999 
liep-ph/9911296 



QED RADIATIVE CORRECTIONS TO e+e // 
WITH REALISTIC CUTS 
AT LEP ENERGIES AND BEYOND 



M. JACK 

DESY Zeuthen, Platanenallee 6, 
D- 15738 Zeuthen, Germany 
E-mail: jack@ifh.de 

After 10 years of steadily increasing the experimental precision at LEP/SLC, there 
is a strong demand on an update of existing programs for fermion pair production. 
We present a rederivation of the 0{a) Bremsstrahlung corrections to e+e~ — > // 
for the semi-analytic program ZFITTER. We focus on observables like total cross 
section and forward-backward asymmetry in the leptonic case with combined cuts 
on acoUinearity angle, acceptance angle, and minimal energy of the fermions. The 
outcome of our analysis is a shift of the predictions by ZFITTER at LEP 1 energies 
off-resonance of a few per mil while at the Z resonance numerical changes can be 
neglected. Thus we obtain for cross sections and asymmetries at LEP 1 a level of 
agreement with other programs of better than per mil, like for the kinematically 
simpler s' cut option. A preliminary analysis of ZFITTER, TOPAZO, and other codes at 
LEP 2 energies showing deviations of several per cent with acoUinearity cuts enforce 
a future examination of higher order effects with different cuts. The predictions 
by LEP/SLC data, however, are not affected within the experimental errors. 



1 Introduction 

In the light of high precision measurements to the Standard Model (SM), fermion pair production 
in e'^e~ annihilation still plays a very important role, e.g. for the extraction of information from 
data on the electroweak symmetry breaking sector of the SM. This impressing increase of precision by 
experiment over the last 10 years was nicely illustrated by m : 



Quantity 


LP 89 (233 events) 


LP 99 (18 X 10^ events) 


Mz (GeV) 
Tz (GeV) 
(light) 


91.17±0.18 

-■^■^3-0.30 

3.0 ±0.9 


91.1871 ±0.0021 
2.4944 ± 0.0024 
2.9835 ±0.0083 



Table 1: Examples for the development of high precision measurements at the Z resonance. 



As one of the most copious processes from LEP to Linear Collider (LC) energies fermion pair 
production also forms an important incoherent background to W-pair production, or more generally, to 
4- fermion final states at higher energies QJ^. Another very interesting application of this 'traditional' 
channel is the possibility to extract limits on extensions of the SM in experiments like LEP or SLC 
when comparing with the theoretical predictions (e.g. Z' and W boson searches, contact interaction 
scales, exchange of leptoquarks, preons, or particles in R-parity breaking supersymmetric models Q). 

We want to focus on high precision measurements to the SM and MSSM. The needed accuracy at 
energies below the Z resonance region from QED Bremsstrahlung, forming the bulk of the radiative 
corrections, is relatively small. However, with the high experimental precision now obtained at LEP 1 
and LEP 2 energies, it is necessary to estimate theoretically these corrections around and above 
the Z resonance region. Especially, hard QED corrections which become more and more important 
with growing cm. energies and resummed soft and virtual higher order effects have to be calculated 

"Talk presented at 14th International Workshop on High Energy Physics and Quantum Field Theory (QFTHEP '99), 
Moscow, Russia, May 27 - June 2, 1999; to appear in the Proceedings. 



precisely taking into account realistic experimental cuts. Having this in mind, the accuracy of the 
theoretical prediction for single contributions to observables should be typically a factor of 10 or so 
better than the experimental error. This is moreover true when considering a high luminosity option 
of a future LC running at the Z resonance ('Giga option) Q. Looking at the experimental situation 
(Table P, we can briefly state that the theoretical predictions therefore have to be better than 0.015% 
at the peak, at the order of 0.03% in the resonance region (w Mz ± 3 GeV), and around 0.5% starting 
at LEP 2 and higher energies. 



2 Realistic observables and the ZFITTER concept 

The semi-analytic approach of the ZFITTER code consists of a fast, one-dimensional numerical 
integration of analytical formulae for different observables like cross sections, asymmetries, polariza- 
tions, and angular distributions with the inclusion of different experimentally relevant cut options 
and resummed 0(a) QED corrections plus dominant higher order effects. The number of analytically 
calculable cuts in such an aproach is of course limited, in our case to two angular cuts and to cuts on 
the final state energies and invariant mass squared. Corresponding numerical programs for fermion 
pair production hke ALIBABA 0, BHM [§, KDRALZ, KK §,|lol, or TDPAZO ||ll| are in this respect com- 
plementary to our approach, as they can in principle treat multi-differential observables with nearly 
arbitrary cuts to the final state phase space, but this at the expense of a clear increase in computing 
time. 

In the ZFITTER approach ||l|,|l|,|l|,^5|j|] , we calculate: 

o{s) ^J— a%.s') pis'/s), (1) 

with s' the invariant mass squared of the final state fermion pair. Cross sections and 

asymmetries like (Tt(s) and Apb[s) are treated in an improved Born approximation, convoluting 
effective Born observables, cr''(s) and A''(s), over s' with a flux function p (radiator) containing the 
photonic corrections. The electroweak and QCD corrections are described with effective couplings in 
the effective Born terms; this is valid due to their smallness at LEP 1 energies. |^ With ZFITTER, three 
different cut options are available: (i) no cut ||l^, (ii) cuts on s' and on the scattering angle d of one 
fermion p^,p^, or (iii) cuts on the fermions' acoUinearity angle, Oacoh on their energies, Ef = Ef, and 



on COST? M, 17|. For a detailed description of the phase space with cuts on final state acoUinearity, 
minimal energies, and/or s' please refer to [p^P]. The effective Born cross sections, (t°(s'), may also be 
chosen according to following approaches: (A) Standard Model, (B) Model Independent, (C) Others 
19 ^. For the latest updates and afs-accounts for the ZFITTER code and other programs please 



consult arapl 20 1 



3 Different codes — comparison and problems 

Focussing first on LEP 1 energies and the s' cut branch (see Section ||), the situation of the ZFITTER 
code up to versions v. 5.x (1998) in comparison with the code TOPAZO can be stated as quite satisfactory. 
At Mz ± 3 GeV the agreement is better than lO^"' [|2^,^. And also at LEP 2 energies and higher we 
can meet the demands by experiment with a deviation of the codes of not more than 1 or 2 per mil 
for different cut values, and with a substantial decrease of this difference below 1 per mil in case of a 
sufficiently large s' cut. This situation does not change when an extra cut on the maximal scattering 
angle cos d is applied [0|o) . 

If we introduced an acoUinearity cut instead of the s' cut, we obtained a comparable agreement 
around the Z resonance as long as we only considered initial state Bremsstrahlung. But as soon 
as we included the initial-final state corrections this agreement deteriorated to 0(3 x 10"'^) which 
grows to an unbearably large discrepancy between the two codes of several per cent at larger energies 
(Vs ~ 100 . . . 200 GeV) [Q. An earfier comparison of the ZFITTER code with the ALIBABA code for the 
s-channel part of the Bhabha scattering branch had already shown similar deviations ||2^,|l^ . 

^ At higher energies, ZZ and WW box corrections contribute at the 1 to 2 per cent level at LEP 2 energies so the 
validity of our non-gauge invariant, effective Born approximation will have to be carefully reexamined there. 
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4 QED Corrections with Acollinearity Cut 



The formulae for photonic 0{a) corrections in ZFITTER with s' cut have been analytically and nu- 
merically multiply checked at LEP energies Jl3| , p^ , p"5| j2|p^ . For the acollinearity cut branch, however, 
there has been no independent check until recently and only little literature available on the exact 
0{a) final state corrections to the total cross section and forward-backward asymmetry [|4j and on 
some formulae related to the initial state corrections (and its combined exponentiation with final 
state radiation) for the angular distribution p5| , p6t . So, a recalculation and documentation of the 
acollinearity cut situation was absolutely mandatory with the main focus first at energies around the 
Z resonance. The slightly more involved Bhabha scattering case with extra t channel contributions is 
kept for a later analysis. 

In our semi-analytic approach it is technically feasible to calculate observables with cuts on one 
of the final state fermions' scattering angles, ??, their energies, Ej their acollinearity angle, 9acoh 
and/or their invariant mass squared, s' . The perturbative calculation of the 0{a) hard corrections to 
total (differential) cross sections then consists of an analytical integration over 3, or respectively 2, 
angles of phase space. The corresponding formulae are finally numerically integrated over s' . 

We observe that neglecting the initial and final state masses at the mentioned high energies 
necessitates a separation of the phase space formed by the remaining two angles of integration, cos 
and cos Oacoh into several different regions. Only where necessary, the masses are kept in order to 
regularize the mass singularities from coUinear radiation of Bremsstrahlung photons. This 'slicing' 
of phase space delivers for each region different analytical expressions for the calculated observables 
[ p^ . For the special cases of either full angular acceptance, i.e. no cut on cos 9, or no cut on the 
acollinearity angle Oacoh the number of different expressions can be substantially reduced and very 
compact formulae can be obtained p^ . This allows one to install a numerically fast second branch 
in the Fortran code when demanding fewer cuts, important for quick data-fitting routines needed by 
experiment. 

The Fortran package acol . f contains a complete collection of the analytical formulae for the 
0{a) corrections and is called from ZFITTER v. 6. 04/06 onwards. The angular distribution will be 
available in versions v.6.2x. A complete collection of all analytical expressions is in preparation p8|. 



5 Results 

The main modifications in the new coding are corrected terms in the QED initial state and interference 
radiator parts. For different regions of phase space, these are different, non-logarithmic contributions 
proportional to cos in the symmetric part of the angular distribution da/dcos d and terms of the type 
a + b cos^ ■& in the antisymmetric part with corresponding changes in the integrated results . 
Omitting these terms in earlier versions is justified with the then anticipated experimental precision 
of only 5 x 10~^ at LEP 1, but not with the higher accuracy now at the Z resonance peak. 

The net effect of the corrected initial state, final state, and interference terms is depicted for ar 
and ApB in Fig. |^ for the energy range 30 to 300 GeV (default flag setting): We compared the cross 
section ratios and the absolute differences of the asymmetries of ZFITTER v. 6. 11 [^, containing the 
new results, with v. 5. 20 H^, still with the old coding, for different acollinearity cut values. The net 
corrections are largest where the radiative return to the Z starts to be prevented by the acollinearity 
cut. For Oacoi < 10° or 25° this sets in at roughly ^/s > 100 GeV, or 115 GeV respectively. These 
corrections to the code are at most roughly 0.5% for ar and 1% for ApB and shrink below 1% at 
higher energies. 

The corrections are mainly due to the new initial-final state interference contributions. This is 
illustrated for a maximal acollinearity angle of 9acoi < 10° in Table ^ where the shifts of the initial 
state corrected total cross sections and asymmetries are shown when switching on the interference 
contributions. One can show that the numerical effects stay below 1% at LEP 2 energies with such a 



cut value |17 20|]. At LEP 1 they are much smaller, of the order of few per mil. 

The corrected initial state and final state terms only have minor effects on ar and ApB and 
amount at most to corrections at the order of 0.1% - 0.2% for ApB at LEP 2 energies for different 
cuts. A detailed analysis of all new modifications to the code can be found in |p^. 
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Figure 1: Net ratios of muon pair production cross sections and differences of forward-bacl<ward asymmetries predicted from 
ZFITTER v.6.11 and v. 5. 20 with three different acollinearity cuts: 6»acol < 10°, 25°, 90°; = 1 GeV; 40° <)9 < 140°; all 

corrections included. 



The numerical comparison of the newly updated ZFITTER version v.6.11 [|6| with TOPAZO's latest 
release version v. 4. 4 ||ll[ now delivers for LEP 1 energies the same high level of agreement as for the 
s' cut (Fig. H): At the peak itself we have a deviation of the codes of 0(10"'') or less for ut and Afb, 
with an acceptable increase to 0(3 x 10~^) for ctt and a slightly worse value for Afb of 0(7 x 10^*) 
with a maximal acollinearity angle of 10°. 

For a wider energy range, we compared ZFITTER v.6.11 (1999) with TOPAZO v.4.3 (1998) and v. 4.4 
(1999) [|O|,0||] and ALIBABA v.2 (1990) in Fig. | An older comparison can be found in g. All 
numbers have been produced with the default settings of the programs. 

At LEP 2 energies the deviation of ZFITTER v.6.11 and TOPAZO v. 4. 4 is at the order of 1% or 
less for different acollinearity cuts and an acceptance cut of 40° < d < 140°. A cross check with 
comparable s' cuts delivers a much better agreement of about 1 per mil or better [p^ . 

In both cases, however, there is a clear peak of the cross section ratios at energies where the 
Z radiative return is not prevented by the cuts. While for the s' cut this discrepancy stays mod- 
erate at the per cent level, it grows up to several per cent for the acollinearity cut. |^ This effect is 
now similar to the observed peak in the ALIBABA - ZFITTER comparison at energies above roughly 
100 GeV. Interestingly enough, when switching off two-loop contributions present in ALIBABA, but 
not in ZFITTER, the agreement improved considerably. Corrections by initial state pair production 
or different exponentiation of initial and final state higher orders, however, do not have a large effect 
here @]. 

Preliminary studies show that a correct description of hard two-loop QED corrections, especially 
for the acollinearity cut option in the ZFITTER code, together with a correct resummation of the soft 
and virtual initial-final state interference contribution, not contained in the ZFITTER code so far, seem 
to play a key role here. Since the acollinearity cut is not as effective in preventing the radiative return 
to the Z as the s' cut, these deviations also survive more profoundly for the acollinearity cut at higher 
energies than for the s' cut. Further comparisons in this respect for both cut options, especially 
including programs KORALZ j|] and KK ||3l[], are intended (see also Q). 



'^The flip of sign of these effects compared to the older versions, TOPAZO v.4.3 and ZFITTER v. 6. 05, is mainly due to a 
corrected interference contribution in the TOPAZO code. Changes to code ZFITTER v. 6. 05 were negligible here. 
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(T^ [nb] with t 


hcol < 10° 








M, - 3 


il/4 - 1.8 




ii4 + 1.8 


il4 +3 




0.21932 


0.46287 


1.44795 


0.67725 


0.39366 


TOPAZO 


0.21776 


0.46083 


1.44785 


0.67894 


0.39491 




7.16 


4.43 


0.07 


+2.49 


+3.17 




0.21928 


0.46284 


1.44780 


0.67721 


0.39360 


ZFITTER 


0.21772 


0.46082 


1.44776 


0.67898 


0.39489 




7.16 


4.40 


0.03 


+2.60 


+3.27 


v4^3 with 0acol < 10° 




M, - 3 


- 1.8 




ii4 + 1.8 


il4 +3 




-0.28450 


-0.16914 


0.00033 


0.11512 


0.16107 


TOPAZO 


-0.28158 


-0.16665 


0.00088 


0.11385 


0.15936 




+2.92 


+2.49 


+0.55 


1.27 


-1.71 




-0.28497 


-0.16936 


0.00024 


0.11496 


0.16083 


ZFITTER 


-0.28222 


-0.16710 


0.00083 


0.11392 


0.15926 




+2.75 


+2.27 


+0.60 


1.03 


1.56 



Tabic 2: A comparison of predictions from ZFITTER v. 6. 11 and TOPAZO v. 4. 4 for muonic cross sections and forward-backward 
asymmetries around the Z peak. First row is without initial-final state interference, second row with, third row the relative 
effect of that interference in per mil. 

6 Summary 

We presented a rederivation of analytical formulae for the 0(a) hard QED Bremsstrahlung corrections 
to e+e~ // in the case of leptonic final states with cuts to the fermions' acollinearity angle and 
energies (/ 7^ e). This was done in the context of the semi-analytic program ZFITTER calculating 
radiatively corrected observables with realistic experimental cuts, e.g. for LEP/SLC applications. 
These corrections to all contributions - initial state, final state, and initial-final state interference - 
are certain non-logarithmic terms which could be neglected earlier but have to be considered now with 
the new high level of experimental precision (e.g. SMz/Mz = 2.2 x 10~^ ||]). 

The older versions of ZFITTER, i.e. versions v. 5. 20 and earlier, derive the 0{a) QED correc- 
tions to ctt with acollinearity cut with a numerical accuracy of about 0.4 % in the Z resonance region 
{Mz ± 3GeV), and similarly for ApB with about 0.13 %. The limiting factor here are mainly the 
initial-final state interference corrections. The new and improved coding in ZFITTER v. 6. 11 Q now 
reproduces the very nice agreement with program TOPAZO v. 4. 4 [|ll| already obtained for the s' cut 
[ll of better than 0.03% (0.1%) for or (Afb) at LEP 1 and better than 0.01% at the Z resonance 
peak. 1^ At higher energies, y/s > 200 GeV, the agreement with s' cut is better than per mil, but 
still only 1 to 2% for the acollinearity cut. Especially in the intermediate energy range, where the Z 
radiative return events are not prevented, this discrepancy peaks and amounts up to several per cent. 
A similar effect is also visible for the s' cut, although to a much lesser extent (below 1%). 

Preliminary studies indicate that higher order corrections, in particular, hard two-loop corrections 
together with an acollinearity cut and a correct resumming of soft and virtual interference effects not 
included in the ZFITTER code seem to be the major underlying cause for this last point | ]3l| , ^ . 
Further studies hereof, also for the angular distribution, and a similar update of the code for the 
Bhabha scattering case are in preparation. 
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